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ABSTRACT 

Normal modes of oscillation of the Sun are useful probes of the solar interior. In this work, we use 
the even-order splitting coefficients to study the evolution of magnetic fields in the convection zone 
over solar cycle 23, assuming that the frequency splitting is only due to rotation and a large scale 
magnetic field. We find that the data are best fit by a combination of a poloidal field and a double- 
peaked near-surface toroidal field. The toroidal fields are centered at r — 0.999i? Q and r — 0.996-R© 
and are confined to the near-surface layers. The poloidal field is a dipole field. The peak strength 
of the poloidal field is 124 ± 17 G. The toroidal field peaks at 380 ± 30 G and 1.4 ± 0.2 kG for the 
shallower and deeper fields respectively. The field strengths are highly correlated with surface activity. 
The toroidal field strength shows a hysteresis- like effect when compared to the global 10.7 cm radio 
flux. The poloidal field strength shows evidence of saturation at high activity. 
Subject headings: Sun: activity, Sun: helioseismology, Sun: magnetic fields 
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1. INTRODUCTION 

Understanding the nature of the Sun's magnetic fields 
— their structure and variability, their generation mech- 
anisms, and their effects on the heliosphere — is one of 
the key aims of current research in solar physics. It is 
generally believed that the magnetic fields are generated 
by a cyclic dynamo that operates somewhere in the solar 
interior. In this paper, we use helioseismology to study 
the global scale internal magnetic fields over the course 
of solar cycle 23. 

Helioseismology is the most powerful tool available to 
solar physicists to study the interior of the Sun. The os- 
cillation frequencies have been used to study the struc- 
ture and dynamics of the solar interior with great preci- 
sion. Magnetic fields, however, have proved to be much 
more challenging. There are a number of important dif- 
ficulties in dealing with magnetic fields in a helioseis- 
mic context. The magnitudes of the signatures in the 
data are quite small, making statistically significant mea- 
surements challenging. Secondly, the interpretation of 
data is very difficult. The physics of wave propagation 
in the presence of magnetic fields is far more complex 
than in the non-magnetic case. Further, the geometry 
of the underlying field strongly affects the signatures in 
helioseismic global mode frequencies, meaning different 
field configurations and strengths can be difficult to dis- 
tinguish from thei r helio seismic signatures. Even worse, 
iZweibel fc Gou eh (1995) showed that because magnetic 
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fields act on mode frequencies both by perturbing the 
thermal structure of the Sun and by changing the wave 
propagation speeds directly, there is a degeneracy be- 
tween magnetic field effects and other thermal perturba- 
tions which cannot be distinguished a priori from helio- 
seismic data. 

Although helioseismic determinations of magnetic 
fiel ds are d ifficul t, there have been many attempts to do 
so. Ilsaald ()1982f ) suggested that the then observed fre- 
quency splittings in the solar acoustic spectrum could 
be cause d by a large scale magnetic field situated in 
the core. iDziembowski fc Goodd (|1984[ ) used an asymp- 
totic approximation to study the effects of magnetic field 
on the splitting coefficients, and IDziembowski &: Goodel 
(1988) argued that a 1 MegaGauss field at the base of the 
convection zone was necessary to explain t he observed 
splitti ng coefficients. However, [Basu (19971); lAntia et all 
(|2000l) placed a limit of 0.3 MG on the field at the base 
of the convection zone; thus the situation was unclear. A 
mega- Gauss magnetic field is also inconsistent with dy- 
namo the^rie^_aiid_TOnstraints from other observations 
C e.g-lD'Silva fc C houdhurilll99l . 

iGough &; Thompson! (|1990f ) developed a formalism to 
compute the effects of rotation and axisymmetric mag- 
netic fields on the freque ncy splittings (disc ussed in the 
following section) , which I Antia et all (|2000f ) used to an- 
alyze the first year of Michelson Doppler Imager (MDI) 
data. They placed limits on the strengths of internal 
toroidal fields, finding a limit of 20 kG at a depth of 
30 Mm, and a limit of 300 kG at the base of the con- 
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vection zone (r = 0.713i? Q ). iDziembowski et all ([2000D 
inverted the mean frequencies and splitting coefficients 
for changes in temperature, and found that the result- 
ing temperature perturbation could be explained by a 
change in magnetic field of 60 kG at a depth of 45 Mm 
( r~0.93.Rf7,). 

IDzie mbowsk i et all ([2001) found that changes in /- 
mode frequencies from solar minimum to solar maximum 
implied a decrease in solar radius with activity, which 
they associated with a change between 4 and 8 Mm in 
depth. In explaining this result with changing magnetic 
fields, they assumed a tangled field, but even so the mag- 
nitude of the change in field strength was strongly depen- 
dent on the radial distribution of the field. The change 
they required was 7 kG for a uniform field, or substan- 
tially le^s_J4_JcG_at_8_^4^n) for an inwa rdly increasing 
field. IChou fc Serebrvanskivl (l200l[2005h looked for sig- 
natures of a change at the base of the convection zone 
from low activity to high activity, and found signs of 
a small change, which they proposed c ould be due to a 
change in mag netic field of 170 - 290 kG. lBaldner fc Basul 
(2008), working with an entire solar cycle's worth of he- 
lioseismic data, found a change in sound speed between 
solar maximum and solar minimum at the base of the 
convection zone, which, if due to a change in magnetic 
field, could indicate a change in field strength of 290 kG 
at that depth. 

In this work, we exploit the fact that we have much 
more helioseismic data than previous investigators had 
access to, and try to get a coherent picture of sub-surface 
solar magnetic field s and their tempor al evolution. We 
extend the work of lAntia et al.l (|2000f ). who considered 
toroidal magnetic fields, to include poloidal fields. This 
means that we can, in principle, consider any axisym- 
metric magnetic field configuration. We compute the ef- 
fects of a wide variety of magnetic field configurations on 
the <Z2 splitting coefficients, and compare them to a so- 
lar cycle's worth of MDI data. It is not clear if the solar 
magnetic field has large scale structure of the form we as- 
sume or whether it is in tangled state due to turbulence 
in the convection zone. Since the effect of magnetic field 
manifests through a quadratic term in magnetic field, 
our estimate may also be applicable to tangled field with 
some degree of approximation. 

2. PERTURBATIONS TO SOLAR OSCILLATION 
FREQUENCIES 

The frequencies of normal modes of oscillation v n £ m 
are degenerate in m in the case of a spherically symmet- 
ric star. Departures from spherical symmetry lift this 
degeneracy. When the departures from spherical sym- 
metry are small, as they are in the case of the Sun, the 
differences in frequency for different values of m will be 
small, and it is natural therefore to express the normal 
mode frequencies in terms of the mean frequency of the 
multiplet v n £ and splitting coefficients af 
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As is common in the current literature, the polynomi- 
als Vj (m) are the Ritzwoller-Lavely formulation of the 
Clebsch-Gordan expansion (jRitzwoller fc Lavelvl Il99l . 



The odd-order splitting coefficients are caused by the ro- 
tation of the Sun, and will not be directly considered 
in this work. The even-order coefficients are caused by 
second order effects of rotation, and by the effects of mag- 
netic fields or any other departure from spherical symme- 
try in the solar structure. In this work, we treat rotation 
and magnetic fields as perturbations on the spherically 
symmetric case, which allows us to avoid explicitly con- 
structing a model of a rota ting, magnetized star . The 
for malism was develope d bv lGough fc Thompson! (|I990D 
and lAntia et al.l (|2000f ) extended the formalism to in- 
clude the perturbation to the gravitational potential (i.e., 
to relax the Cowling approximation) and to include dif- 
ferential rotation. 

The first order correction to the mode frequencies due 
to rotation affects only the odd-order splitting coeffi- 
cients. These effects are due to the perturbation of the 
mode frequencies by advection of the waves. The second 
order correction affects only the even-order splitting co- 
efficients, and is caused by the perturbation to the eigen- 
functions and the centrifugal force. The odd-order coeffi- 
cients can be used to determine the rotati on profile £l(r) 
(|Thompson et al] [l996; Sch ou et al.lll998T ). which can in 
turn be used to compute the second order rotation cor- 
rection (|Antia et al.ll2000f ) to the even-order coefficients. 
This correction needs to be made if the magnetic per- 
turbation is comparable in size to second order rotation 
effect, which appears to be the case (jGough fc Thompson! 
[199(3: lAntia etan[200l . 

In this work, we consider two different axisymmctric 
magne tic field configurations : toroi dal and poloidal. Fol- 
lowing [Gough~2jrhom2ioii| (|1990f ). the toroidal field is 
expressed in the form 



0,0,o(r)— P k (cost 
do 
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where are the Legendre polynomials of degree k and 
o(r) describes the radial profile of the magnetic field. We 
consider only even values of k to ensure antisymmetry 
about the equator, consistent with the observed field at 
the surface. The poloidal field is assumed to be of the 
form 
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where b(r) describes the radial profile of the magnetic 
field. In this case we use only odd values of k to ensure 
that the field is antisymmetric about the equator. With 
appropriate combinations of these two fields we can, in 
principal, represent any axisymmetric magnetic field. 

The effect of these magnetic field configurations on the 
frequency split tings of p-modes is calcula te d using the 
formu lation of IGough fc Thompson! (|1990D ; lAntia etlrtl 
(1200(3) . There are two ways in which the magnetic field 
can affect the frequencies, one is the so-called direct effect 
due to the additional force, and the second is the distor- 
tion effect due to the equilibrium state being distorted 
from the original spherically symmetric one. Both these 
effects are included in all calculations. These formula- 
tions treat the effect of these magnetic fields separately. 
Unfortunately, the effect of magnetic fields is not linear 
and hence strictly the contributions from two different 
configurations cannot be added. In principle, there will 
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be some cross-terms when the combination of toroidal 
and poloidal fields have a region of overlap in the solar 
interior. In this work, we neglect these terms and add 
the contributions from toroidal and poloidal fields to get 
the total effect. We expect the cross terms to be small. 

3. DATA 

The data we use for comparison are 72-day mode pa- 
rameter sets from the Michelson Doppler Imager (MDI) 
on the SOlar and Heliospheric Observatory (SOHO). We 
use mode p arameter sets from the corrected pipeline de- 
scribed by iLarson fe Schoul (120081) . The original MDI 
analysis pipeline (|Schoulll999f ) did not take in to account 
a number of instrumental effects which introduced secu- 
lar trends in the mode parameter sets. In particular, the 
plate scale of the MDI instrument has changed somewhat 
over SOHO's mission, and this results in an apparent 
change in th e solar radius if n ot pr operly corrected in 
the analysis. [Baldncr fe Basul (|2008T ) found a signature 
in the mean frequencies which became increasingly sig- 
nificant over the course of the solar cycle. A repetition of 
that work with reanal yzed mode paramete r sets removed 
this effect completely (|Baldner et al.l l2008h The splitting 
coefficients, which we focus on in this work, suffer from 
similar instrumental effects as the mean frequencies, and 
hence we use the reanalyzed data in this work. 

We include 56 mode sets which cover solar cycle 23. 
The mode sets arc identified by the MDI start day, begin- 
ning with set #1216 (start day 1 May 1996), and ending 
with set #5320 (start day 27 July 2007). The coverage 
begins and ends at low activity, with a 10.7 cm radio flux 
of 72.7 SFU for the first set and a flux of 69.1 SFU for 
the last set. The highest activity set, #3160 (start day 
27 August 2001), has a 10.7 cm flux of 223.9 SFU. 

We fit only the splitting coefficients, as the higher 
order splitting coefficients have larger errors, and as such 
did not distinguish well between different field configu- 
rations. The rotation profi le determined from the odd- 
order splitting coefficients (|Antia et alj l2008j) was used 
to calculate the second-order contribution to the even- 
order coefficients, and this contribution was subtracted 
from the data. 

4. RESULTS 
4.1. Models 

In Fig. [TJ we show the second splitting coefficient for 
four different poloidal field configurations. The actual 
quantities plotted are la-x-, both as a function of frequency 
v and as a function of the lower turning radius of the 
modes, r t . The radial profile in this case is taken to be 

b(r) = B Q r- k , (4) 

where Bq is a constant which determines the peak field 
strength and r is the radial distance measured in units of 
solar radius. The models shown in Fig. [T] all have a peak 
strength of B = 1 G (note that they do not all have the 
same value of Bq). The most obvious difference between 
different order poloidal fields is that for the fc — 1 field 
the splitting coefficients are all positive, whereas for the 
higher order fields they are largely negative, although the 
shallow modes have positive a,2- 
The toroidal field we employ is similar to that used by 
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Fig. 1. — Splitting coefficients lai due to poloidal magnetic fields. 
The left hand panels are shown as a function of frequency u, the 
right hand panels are shown as a function of lower turning radius rt ■ 
The four configurations shown have peak field strengths of 1 G at 
the surface. The fields have four different values of k. To facilitate 
direct comparison with later figures, only modes measured in the 
MDI data (specifically, the high activity set #3160) are plotted. 

lAntia et al.l (|2000f ). with a radial profile given by 

a(r) = J V8^(l~(^) 2 ) if>-r |<<*o (5) 
[ otherwise 

where p is the gas pressure, /3o is the ratio of the mag- 
netic to gas pressure at ro, and ro and do are position 
and width of the field. As is the case for the poloidal 
fields, the toroidal field corrections are linear in magnetic 
field strength squared. Excepting field strength, there- 
fore, our toroidal fields are described by three quantities: 
the order of the Legendre polynomial fc, which deter- 
mines the latitudinal distribution of the field, the cen- 
tral radius ro, which determines the location, and the 
width do. Figure [2] shows the splitting coefficients due 
to toroidal fields with different values of fc but the same 
radial profile (in this case, (3q = 10~ 4 , r = 0.999i? Q and 
do = O.OOli? ). For the 02 coefficient, the order k of the 
field makes very little difference except to the scale of 
the perturbation — increasing fc for the same [3q effec- 
tively increases the total amount of flux, but except for 
this effect, the 0,2 coefficients are not sensitive to different 
latitudinal distributions. For the remainder of the work, 
therefore, we restrict ourselves to fc = 2 fields. 

Figure [3] shows the splitting coefficients £02 for near- 
surface toroidal fields with different central radii ro and 
widths do as a function of frequency. Figure [4] shows 
the same, but as a function of the lower turning ra- 
dius, r t . The behavior of the splitting coefficients is 
not surprising. In general, the fields which penetrate 
below the surface show oscillatory behavior as a func- 
tion of freque n cy similar to that seen in m ode frequency 
(iGoughl fl990L iGoueh fc Thompson! [199(1 a nd used by 
iRoxburgh fc Vorontsovl (|1994fk iBasu et all (|1994fl and 
others to study the convection zone base. The period of 
these oscillations is related to the acoustic depth of the 
perturbation in the structure. Decreasing the depth of 
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Fig. 2. — Splitting coefficients tai due to toroidal field with different latitudinal distributions. The upper panels show the coefficients as 
a function of frequency u, the lower panels show the coefficients as a function of lower turning radius rt- All the results are with fto = 10~ 4 , 
ro = 0.999-Rq and dg = 0.001-Rq. Only modes present in the MDI data have been plotted. 



the perturbation lengthens the period of the oscillatory 
behavior. Fields which are confined near the surface, on 
the other hand, do not exhibit oscillatory behavior, but 
instead resemble the 'surface term' correction which is 
removed in structure inver sions (e.g.. fDziembowski et all 
I199CA lAntia fc Basul 119941 ). Increasing the width of the 
perturbation smears out the oscillatory signature, as seen 
in Fig. [3] Because all the modes sampled have lower turn- 
ing radii below the magnetic fields considered here, there 
are no obvious signatures in the splitting coefficients as 
a function of r t . 

In addition to fields near the surface, in Fig. we 
show the splitting coefficients due to some toroidal fields 
located at the base of the convection zone. The fields 
shown differ only in the width do of the fields. Unlike 
the surface fields shown in previous figures, the deep field 
signatures show both positive and negative splitting coef- 
ficients. These models are most interesting as a function 
of lower turning radius r t . The splitting coefficients are 
positive above the center of the magnetic field, and neg- 
ative below the center of the magnetic field. Further, 
as the width of the field is increased, the width of the 
perturbations to the splitting coefficients (in r t figure) 
increases as well. 

The 04 splitting coefficients due to various poloidal and 
toroidal fields are shown in Fig. [51 which shows the results 
for two poloidal fields, the k = 3 and k — 7 fields, as well 
as two toroidal fields with different values of k (k = 2 and 
k = 8), each with r = 0.999i? Q and d = 0.001i? Q . The 
k = 1 field has essentially no effect on the 04 splitting- 
coefficients. 

4.2. Fits to observed data 

In order to choose the fields which best match the ac- 
tual data, we have computed the splitting coefficients for 
a large grid of field configurations, with fields through- 
out the convection zone. For poloidal fields we varied k 
— the form of the radial profile was found not to mat- 
ter very much for the splitting coefficients, so long as the 
field penetrated below the surface. For the toroidal fields, 



we varied the location ro, the width of the field do, and 
the latitudinal distribution with k. The range in ro was 
between 0.70-R© and I.ORq. The values for do ranged 
from 10~ 4 i? Q to 0.2i? Q . In order to judge goodness-of- 
fit, we use the x 2 statistic. For both the poloidal and 
the toroidal fields, the perturbations vary linearly with 
the square of the field strength, so to fit the field, we 
allowed the field strength to vary freely, and chose the 
strength that minimized the x 2 ■ We have computed the 
X 2 for all the field configurations in our grid, as well as 
for many combinations of two and three different fields. 
The results we present below represent the best fits from 
the entire grid of computed models. 

The largest signal-to-noise ratio in a-i is found at peak 
activity, and so the highest activity set ought to be the 
easiest to fit. Comparison of different field configura- 
tions with the splitting coefficients at high activity are 
shown in Fig. [3 and the fits are shown both as a function 
of frequency and as a function of lower turning radius. 
The residuals, normalized by the errors in the data, are 
also shown. A fit to a k — 1 poloidal field is shown 
in panel (a). The reduced x 2 f° r this fit is 16, and it 
is evident that the field does a poor job of reproducing 
the observed splitting coefficients. Higher order poloidal 
fields are considerably worse, as an examination of Fig. [T] 
will show — these fields perturb all the splitting coeffi- 
cients negatively, whereas the observed splittings are all 
positive. Panel (b) shows the effect of a toroidal field 
situated near the surface. Although we attempted to fit 
toroidal fields throughout the convection zone, fields not 
located very near the surface were extremely poor fits 
to the data. The field shown in panel (b) is the best 
fit for a single toroidal field, with ro = 0.999-R© and 
do = O.OOli?©. The reduced x 2 is 5. The residuals are 
mostly without structure in r t , but are oscillatory in fre- 
quency, a hint that there could be a second, somewhat 
deeper field. The splitting coefficients at peak solar ac- 
tivity cannot be well fit by either a toroidal field or a 
poloidal field of the form considered by us. 

The third field configuration shown (panel c) is a com- 
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Fig. 3. — Splitting coefficients l<i2 due to near-surface toroidal magnetic fields, as a function of frequency v. The results are shown 
for fc = 2 with five different values of central radius ro (from 0.996-Rq to Rq), and three different values of the width of the field do 
(0.0001, 0.001, and 0.005)R Q . Only modes present in the MDI data have been plotted. 




Fig. 4.— 



Same as Fig. [3] but plotted as a function of the lower turning radius rt. 
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Fig. 5. — Effects of toroidal magnetic fields at the base of the con- 
vection zone on the la,2 splitting coefficients. Results for two mag- 
netic field configuration with k = 2, ro = 0.71-Rq and (3q = 10~ 4 
are shown. The fields have widths of do = O.OliJ© and do = 0.1-Rq 
for the top and bottom panels, respectively. Left hand panels show 
the splitting coefficients as a function of frequency v, right hand 
panels show the splitting coefficients as a function of lower turning 
radius rt- Only observed modes have been plotted. 
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Fig. 6. — Effects of various magnetic fields on the £04 splitting 
coefficients, as a function of both frequency v (left hand panels) and 
lower turning radius rt (right hand panels). The top two panels 
show the results for poloidal fields with k = 3 (panel a) and k = 7 
(panel b). The bottom two panels show the results for toroidal 
fields, both with /3 = 10~ 4 , r = 0.999-Rq and do = O.OO1H . 
Panel (c) is for a k = 2 field and panel (d) for a k = 8 field. Only 
observed modes have been plotted. 

bination of a k — 1 poloidal field and a near-surface 
toroidal field (r = 0.999i? Q , d = O.OOli?© — the same 
field from panel (b)). This combination of fields yields 
a much better fit to this data set, with a reduced \ 2 of 
2.8. Using a surface toroidal field instead of a poloidal 
field does not fit the data as well, although it is an im- 
provement over the fit in (b), with a x 2 value of 3.5. 
Like the toroidal-only fit, the residuals are more or less 
without structure in r t , but show oscillatory behavior in 
frequency. The peak field strengths of the two fields are 
133 G and 368 G for the poloidal and toroidal fields, re- 
spectively. The residuals from this fit can be fit by a 
toroidal field centered at r — 0.996i? Q , so in panel (d) 
we show the best fit to the data: a k = 1 poloidal field 
with two toroidal fields, one centered at r = 0.999i? Q 



and another centered at ro = 0.996i?Q. Both toroidal 
fields are k = 2 fields and have widths do — 0.001i? Q . 
The poloidal field has a peak field strength at the sur- 
face of 124 ± 18 G, while the toroidal fields have peak 
field strengths of 380 ± 30 G and 1.4 ± 0.2 kG, respec- 
tively. The reduced x 2 °f this fit is 1.7. Attempts to 
fit the data with a single toroidal field which occupies 
the same region as the two field in this fit did not yield 
a good fit — the data seem to require a double peaked 
field. 

Figure [8] shows the 04 splitting coefficients for two dif- 
ferent data sets (#2224 and #3160), and the same mod- 
els shown in Fig. [7] panel (d). The errors on the data 
are large compared to the signal — thus the normalized 
residuals are comparable to those for the 0,2 coefficients, 
but some other field configurations also fit the 04 coeffi- 
cients equally well, so we do not use them to constrain 
the field configurations or determine the field strengths. 

Having fit the high activity set, we repeat the fits for all 
56 sets in our study. Figure[9]shows fits of the same k — 1 
poloidal plus toroidal field combination to six representa- 
tive mode sets, covering the rise and fall of solar cycle 23. 
The first set, #1216, is the first 72-day mode set from the 
MDI program, and is a low activity set. Two rising phase 
sets are shown, #2224 and #2728, with 10.7 cm fluxes 
of 131.4 SFU and 187.3 SFU, respectively. The high ac- 
tivity set from Fig. [7] is shown, and a declining phase 
set (#3952, -Fio.7 = 126.4) and a set from the current 
minimum (#4744). For all sets, the same combination 
of poloidal and toroidal fields — but with differing field 
strengths — was found to be the best fit. The residuals 
show the same structure with frequency. 

No mode set in our study is well fit by any magnetic 
field at the base of the convection zone. For low activity 
sets (10.7 cm flux of less than 100 SFU), we can find an 
upper bound on the field that could be present in the 
data. For a field centered at the base of the convection 
zone with k = 2 and do = 0.01i?o, we find that a fields 
of up to 300 kG can be fit to the data, although this is 
an upper limit, not a detection, since the models of that 
field strength or lower give the same \ 2 as a zero field 
strength model. At high activity, the dominant signal 
is from the surface, which we have attempted to explain 
with magnetic fields located in those layers. 

The field strengths of the poloidal and shallow toroidal 
field fits to all 56 sets used in this study are shown in 
Fig. [10] Also shown in this figure is the ratio of the 
poloidal field strength to the ro = 0.999i?o toroidal field 
strength. With the exception of the low activity sets 
at the beginning and end of the solar cycle, where the 
uncertainty in the fits is relatively large, the ratio be- 
tween the poloidal and toroidal field strengths is roughly 
constant. The field strengths from Fig. \W\ are corre- 
lated with global activity indices from solar cycle 23. 
The correlation coefficients are 0.90, 0.93, and 0.92 for 
the poloidal and two toroidal field components, respec- 
tively. In Fig. [TTJ we plot the toroidal and poloidal field 
strengths as a function of one such global index, the 
10.7 cm radio flux. The field strengths prove to be highly 
correlated with activity, although there is a hysteresis- 
like effect evident in the toroidal field strengths — the 
rising phase (shown in blue) is weaker than the declin- 
ing phase fields. The same effect may also be present at 
low activity in the poloidal field strengths. The poloidal 
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Fig. 7. — Fits to observed splitting coefficients lai for different magnetic field configurations. Four different fits are shown, both as a 
function of frequency v and as a function of lower turning radius rj. The data are shown in black, and the modeled points in red. The 
residuals, scaled by the errors in the data, are also shown below the comparisons. The data are the la,2 splitting coefficients from an MD1 
72 day mode parameter set, taken at the peak of solar cycle 23 (MDI set #3f 60, start day 200f Aug 27). Panel (a) shows a fit from a k = 1 
poloidal field. Panel (b) shows the fit from a near-surface toroidal field. Panel (c) shows the best fit field combination with two fields to 
this data set — a combination of a dipole poloidal field with a toroidal field located just below the surface (rg = 0.999-Rq, do = O.OOI.R0). 
Panel (d) shows the best fit with three fields — the same poloidal and toroidal field as in panel (c) (though with slightly different field 
strengths) and another toroidal field at t*o = 0.996-Rq and a width of do = O.OOfi?0. 



field strengths do seem to saturate at high activity. The 
strengths of the two toroidal fields are extremely well 
correlated. 

5. DISCUSSIONS AND CONCLUSIONS 

We have attempted to use the first even order splitting 
coefficient (02) to infer the configuration and strength of 
the Sun's internal magnetic fields over the course of solar 
cycle 23, assuming that the entire signature in 0,2 after 
correction for rotation effects is magnetic and that the 
fields are axisymmetric. The field that we have found is 
a combination of poloidal field and a double-peaked near- 
surface toroidal field. The strengths of the poloidal and 
toroidal components, at least for high activity period, are 
well correlated. The relative strengths of the two toroidal 
fields are also extremely well correlated. 

Although the fits we have shown are the best fit to the 
data from the grid of models that we have computed, we 
can say nothing about the uniqueness of these fits over 
the set of all possible magnetic field configurations in the 
solar interior. In particular, the choice of radial profile of 
the toroidal fields is virtually limitless, and by restrict- 



ing our work to profiles of the form ([2]) , we have limited 
our search to a restricted class of fields. It is possible 
that there are fields we did not consider with quite dif- 
ferent radial and latitudinal distributions which fit the 
data as well as the fields we have presented as best fits. 
In addition, as noted above, it is not strictly correct to 
add the splitting coefficient perturbations together as we 
have done without explicitly accounting for the pertur- 
bations arising from the cross terms. We do not expect, 
however, these corrections to be significant, and a full 
treatment of these corrections would be considered in a 
future work. 

Our inferred magnetic field does not change its latitu- 
dinal distribution over the course of the solar cycle. This 
is in part due to the fact that we are only fitting the 
(Z2 coefficient (as noted above, the higher order splittings 
had large errors), so our sampling of the interior is not 
really latitudinally sensitive. Thus, we do not see a but- 
terfly diagram in our magnetic fields. lUlrich fc Bovdenl 
(2005) measured the surface toroidal component of the 
solar magnetic field over almost an entire 22-year cycle. 
The field they measure is roughly a tenth of the peak 
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Fig. 9. — Fits to measured splitting coefficients la,2 for six different sets throughout solar cycle 23. The data (shown in black) are from 
MDI 72 day mode parameter sets. The magnetic field configuration is the same as panel d in Figure [7] a dipole poloidal field and two 
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Fig. 10. — The strength of the inferred magnetic fields as a func- 
tion of time over solar cycle 23. Each MDI 72 day mode parameter 
set is fitted by with the same magnetic field configuration as Fig. [9] 
The strengths (in Gauss) of the poloidal field at r = 0.999i?Q (solid 
black line) and the toroidal field at r = 0.999-Rq (dashed line) are 
shown in the upper panel (a). The middle panel (b) shows the 
same quantities as in the upper panel, but this time at a radius of 
r = 0.996Rq. The lower panel (c) shows the ratio of the poloidal 
field strength to the toroidal field strength at r = 0.999Rq. The 
ratio of poloidal to toroidal at r = 0.996Rq looks very similar. 




Fig. 11. — The strength of the inferred magnetic fields as a func- 
tion of 10.7 cm radio flux. The top panel shows the ro = 0.999-Rq 
toroidal field strength, the middle panel shows the r = 0.996-Rq 
toroidal field strength, and the bottom panel shows the poloidal 
field strength at r = 0.999-Rq. Rising and declining phase are dis- 
tinguished with blue circles for the rising phase and red triangles 
for the declining phase. The toroidal field shows a hysteresis effect. 
The poloidal field shows some hysteresis at low activity, as well as 
a hint of saturation at high activity. 

strength of our toroidal field. Strictly speaking, how- 
ever, we see no toroidal field at all at the surface, since 
in our inferred field, the field strength becomes zero pre- 
cisely at r = IRq. The peak strength that we measure, 
however, is only 700 km below the sur face, and the field 
could penetrate the surface somewhat. lUlrich fc Bovdenl 
( 2005) find a field which gives a /3 ~ 6 x 10 5 at the sur- 
face at high activity, and drops to nothing at low activity, 
while we find a field that changes from /3q ~ 10~ 4 at low 



activity to /3o — 2 x 10~ 3 at high activity at a radius of 
r = 0.999-R© (a depth of approximately 700 km). 

Recently, attention has been focused on the strength 
and configuration of the quiet Sun surface magnetic field. 
IHarvev et alj ([2007) reported the presence of a 'seething' 
horizontal magnetic field with an rms field strength of 
1.7 G. With the launch of Hinode (Solar B), the high 
spatial resolution of the onboard spectropolarimcter has 
been used to study the horizontal fi elds of the solar pho- 
tosphere. iLites et all (|2007l I2008f ) have measured the 
horizontal flux, which they find to be 55 G, compared 
to the average vertical flux of 11 G. IPetrie fc Patrikeeval 
(2009) found that the zonal component (component in 
the East- West direction) was much smaller than the ra- 
dial component, reporting an inclination angle of less 
than 12° from vertical in the East- West direction. 

The fields being studied in the aforementioned works 
are generally very tangled fields which thread through 
the intergranular lanes and so they are not axisymmet- 
ric fields. It is worthwhile to compare our results with 
theirs, since tangled fields on local scales can organize 
into roughly axisymmetric fields on global scales. How- 
ever, the contribution to splittings are more sensitive to 
(B 2 ) rather than (B) 2 and hence tangled field may also 
contribute to it, even when the average (B) is very small. 
Further, considering the general behavior of the pertur- 
bation to the mode frequencies, our inference about the 
location of required magnetic field is more robust as a 
different location will yield a very different behavior of 
splitting coefficients. The exact magnitude of the field 
may depend on the assumption of geometry and on it 
being tangled or large scale. Nevertheless, we believe 
that our estimate is of the right order, though the statis- 
tical errorbars obtained by us may not be realistic. The 
systematic errors in these estimates would be certainly 
larger. The domin ance of poloidal field orien tation at 
the surface found by IPetrie fc Patrikeeval (|2009f ) is found 
in our own results — at the surface, the toroidal field 
is weak or vanishing, b ut the polo i dal field rem ains. In 
the period analyzed bv ILites etaTl (|2007L l2008h . we find 
a poloidal field strength of 40 G, and a toroidal field of 
90 G at a depth of 700 km. The vertical flux they find 
(11 G) is weaker than what we detect, but their 55 G hor- 
izontal flux may be roughly consistent with our toroidal 
field. 

ISchrTi ver fc L iu (200cf) found that the dipole moment 
of the surface magnetic field, measured from MDI mag- 
netograms, was half the strength in 2008 that it was in 
1997, during the last solar minimum. We do not see such 
a difference from the beginning of our period to the end 
in fact, we find the poloidal field strength is slightly 
higher during the current minimum, although the level 
of the difference is within the errors, and our data sets 
end in 2007, so the comparison is not contemporaneous. 

Hysteresis in the relations between activity indices 
has been observed before ; for example in the rela- 
tion between low d e gree flAnguera Gubau et all IT992T: 
iJimenez-Reves et al.l Il998l) and intermediate degree 
( Trip athy et al.l 120 00. 2001) acoustic modes and global 
magnetic indices. It should be noted that an analysis of 
a full solar cycle's worth of intermediate degree p-modes 
data does not show any hyster esis in mean frequencie s 
as a function of 10.7 cm flux (jBaldner fc Basull200l . 
iTripathv et al.l (|2001f ) noted that, among the global mode 
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indices, the relation between global line-of-sight mag- 
netic flux and 10.7 cm radio flux showed a hysteresis 
effect, but the relation between the radiative indices and 
10.7 cm flux did not. iMoreno-Insertis fc Solankil (|2000l ) 
argued that the observed hysteresis could be almost en- 
tirely due to the latitudinal distribution of magnetic flux 
on the surface of the Sun. We believe that this is a 
compelling explanation for the hysteresis that we find. 
The 10.7 cm flux is the integrated flux received at the 
Earth and does not contain any information about the 
latitudinal variation, while the ci2 splitting coefficient is 
associated with definite lat itudinal variation, given by 
P2(cos#) (|Antia et al. 1 12001). and hence the two would 
not be the same. More importantly, we expect the ac- 
tual magnetic fields in the near surface layers to drift 
equatorward — as the surface fields do. 

Few conclusions can really be drawn from this work 
with respect to dynamo theory since the fields we have 
inferred are predominantly shallow fields, whereas most 
dynamo models operate much deeper down, in the shear 
layer at and below the base of th e convection zone . 
(some useful recent reviews include |Ossendriiverl 120031 : 
iCharbonneaul 120051 : iMiesch fc Toomrd l2009f ). The up- 
per limits that we place on fields at that dept h are 
consi s tent with earlier helioseismic results (e.g. , iBasul 
119971 lAntia et all 120001: iChou & Serebrvanskivl 120021 . 



120051 iBaldner fc Basul I2008Q . Many deep-seated dy- 
namo mechanisms predict an anticorrelation between 
the poloidal and toroidal field components, as the dy- 
namo converts poloidal to toroidal field and toroidal field 
back to poloidal. We do not see any evidence of such 
conversion. Some dynamo mechanisms, however, op- 
erate in the near-surface shear layer (e.g. iBrandenburd 
2005). Although the fields generated in these models 
are generally extremely tangled, on global scales these 
fields can have to r oidal and poloidal components (e.g. 
IBrown et all 120071 1200H) . In particular, although they 
were considering a more rapidly rotating star than the 
Sun, IBrown et al.l (|2007l ) noted that their field contained 
both a poloidal and a toroidal component, and that the 
toroidal component was much the stronger of the two. 



This work utilizes data from the Solar Oscillations In- 
vestigation/ Michelson Doppler Imager (SOI/MDI) on 
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